B AND D MESONS IN LATTICE QCD 



FERMILAB-CONF-00/256-T 



ANDREAS S. KRONFELD 
Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia, IL, USA 

Computational and theoretical developments in lattice QCD calculations of B and D mesons are 
surveyed. Several topical examples are given: new ideas for calculating the HQET parameters A and 
Ai; form factors needed to determine |V c j,| and |V^i,|; bag parameters for the mass differences of the 
B mesons; and decay constants. Prospects for removing the quenched approximation are discussed. 



1 Introduction 

In the standard model, interactions involv- 
ing the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix violate CP, with strength propor- 
tional to the area of the "unitarity triangle." 
Fig. ^ shows a recent summaryfcl of the tri- 
angle. The dominant uncertainties are theo- 
retical, coming from non-perturbative QCD. 
Each blob shows experimental uncertainties 
for fixed theoretical inputs, so the range of 
blobs illustrates the theoretical uncertainties. 
If measuring the apex [p, fj) were the only 
goal, one might conclude from Fig. [I] that the 
most pressing issue is to reduce the theoreti- 
cal uncertainties, which would require greater 
investment in computing for lattice QCD. 

Measuring (p, fj) is not the most excit- 
ing goal, however. "High-energy physics is 
exciting and will remain exciting, precisely 
because it exists in a state of permanent 
revolution. "□ That means we would prefer to 
discover additional, non-KM sources of CP 
violation. Indeed, "it is possible, likely, un- 




P 



Figure 1. Constraints on the unitarity triangle from 
CP conserving B decays and indirect CP H°l a t,i° n 
in the kaon. From Plaszczynski and Schune.H 



avoidable, that the standard model's picture 
of CP violation is incomplete. "13 

Lattice QCD can aid the discovery of 
new sources of CP violation and may be es- 
sential. A lot of information will be neces- 
sary to figure out what is going on at short 
distances. One way to think about this is 
sketched in Fig. ||. The triangle BjA is deter- 
mined from (quark-level) tree processes. The 
side B requires \V u d\ from n — > pe~D, and 
\V ub \ from B~ -> p°l~P or B° -> n+l'D; 
the angle 7 requires the CP asymmetry of 
B ± -> D^pK* (or B s -> DfHJ); the side 
A requires \V c d\ from D° — > n~l + v, and \V c b\ 
from B~ — » D°^l~v. One could call this the 
"tree triangle" . The triangle aC{3 is deter- 
mined from mixing processes (including in- 
terference of decays with and without mix- 
ing). The angle a requires the asymmetry 
of B — > p7r; the side C requires \V t d\ from 
Am B0 , and |Vt&| from t — ► W + b; the angle 
f3 requires the asymmetry of B — > J/ipKs- 
One could call this the "mixing triangle" . 

Checking whether the mixing triangle 
agrees with the tree triangle tests for new 
physics in the amplitude of B^-B® mixing. 
New physics in the magnitude muddles the 
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Figure 2. Two different unitarity triangles: the tree 
triangle B^A, and the mixing triangle aCj3. 
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extraction of \Vtd\, and new physics in the 
phase muddles the extraction of angles a 
and (3. Similarly, taking Am B0 , B s — > 

DfK*, and B s -> J/iprj^ (or J/ip<j>) sorts 
out new physics in B®-B® mixing. 

These tests are impossible without know- 
ing the sides accurately, so hadronic matrix 
elements are needed. In a few cases a symme- 
try provides it, e.g., isospin cleanly yields the 
matrix element for n — > pe~D. For the oth- 
ers, we must "solve" non-perturbative QCD 
and, therefore, you need lattice calculations. 

You are probably tired of waiting and 
may ask why results should come any time 
soon. The fastest of today's computers are 
now powerful enough to eliminate the sorest 
point: the quenched approximation. Also, 
(lattice) theorists have slowly developed a 
culture of estimating systematic uncertain- 
ties, which is now not bad and would im- 
prove if more non-practitioners became suffi- 
ciently informed about the methods to make 
constructive suggestions. 

The rest of this talk starts with some the- 
oretical aspects that might make it easier for 
the outsider to judge the systematic errors of 
heavy quarks on the lattice. Then I show re- 
cent results needed for the sides A, B, and C,. 

2 Lattice Spacing Effects (Theory) 

Lattice QCD calculates matrix elements by 
computing the functional integral, using 
a Monte Carlo with importance sampling. 
Hence, there are statistical errors. This part 
of the method is well understood and, these 
days, rarely leads to controversy. When con- 
flicts do arise, they usually originate in the 
treatment of systematics. The non-expert 
does not need to know how the Monte Carlo 
works, but can develop some intuition of how 
the systematics work. Don't be put off by lat- 
tice jargon: the main tool is familiar to all: 
it is effective field theory. 

Lattice spacing effects can be cataloged 
with Symanzik's local effective Lagrangian 



(LE£)fl Finite-volume effects can be con- 
troled and exploited with a general, mas- 
sive quantum field theory.i The computer al- 
gorithms work better for the strange quark 
than for down or up, but the dependence on 
m q can be understood and controlled via the 
chiral Lagrangian.Q Finally, discretization ef- 
fects of the heavy-quark mass ttiq are treated 
with HQET2 or NRQCD.1 In each case one 
can control the extrapolation of artificial, nu- 
merical data, if one generates numerical data 
close enough to the real world. 

Volume effects are unimportant in what 
follows, and chiral perturbation theory is 
a relatively well-known subject. Therefore, 
here I will focus on the effective field theories 
that help us control discretization effects. 

2.1 Symanzik's LEC 

Symanzik's formalism^ describes the lattice 
theory with continuum QCD: 

£lat = £ CO nt + 'Y^a Si C i {a;n)O i {^), (1) 

i 

where the symbol = means "has the same 
physics as" . The LE£ on the right-hand side 
is defined in, say, the MS scheme at scale 
fi. The coefficients Oj describe short-distance 
physics, so they depend on the lattice spac- 
ing a. The operators do not depend on a. 

If Aqcdo is small enough the higher 
terms can be treated as perturbations. So, 
the a dependence of the proton mass is 

m p (a) = m p + aC a F(p\'ipa- ■ Fip\p), (2) 

taking the leading operator for Wilson 
fermions as an example. To reduce the sec- 
ond term one might try to reduce a greatly, 
but CPU time goes as a~( 5 or 6 ). It is more 
effective to combine several data sets and ex- 
trapolate, with Eq. (||) as a guide. It is even 
better to adjust things so C a F is 0{ar s ) or 
0(a), which is called Symanzik improvement 
of the action. For light hadrons, a combi- 
nation of improvement and extrapolation is 
best, and you should look for both. 
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2.2 HQET for large m Q 

The Symanzik theory, as usually applied, as- 
sumes m q a <C 1. The bottom and charm 
quarks' masses in lattice units are at present 
large: m^a ~ 1-2 and m c a about a third 
of that. It will not be possible to reduce a 
enough to make m^a <C 1 for many, many 
years. So, other methods arc needed to 
control the lattice spacing effects of heavy 
quarks. There are several alternatives: 

1. static approximation^] 

2. lattice NRQCD0 

3. extrapolation from rriQ w m c up to mj, 
3'. combine 3 with 1 

4. normalize systematically to HQET0 

All use HQET in some way. The first two 
discretize continuum HQET; method 1 stops 
at the leading term, and method 2 carries 
the heavy-quark expansion out to the desired 
order. Methods 3 and 3' keep the heavy 
quark mass artificially small and appeal to 
the \jvfiQ expansion to extrapolate back up 
to mj. Method 4 uses the same lattice action 
as method 3, but uses the heavy-quark ex- 
pansion to normalize and improve it. Meth- 
ods 2 and 4 are able to calculate matrix ele- 
ments directly at the 6-quark mass. 

The methods can be compared and con- 
trastediry describing the lattice theories with 
HQET.E3 This is, in a sense, the opposite of 
discretizing HQET. One writes down a (con- 
tinuum) effective Lagrangian 

Aat = C lat ( m Q a i M)0HQ ET fci), ( 3 ) 
n 

with the operators defined exactly as in con- 
tinuum HQET, so they do not depend on 
toq or a. As long as rriQ Aqcd this de- 
scription makes sense. There are two short 
distances, 1/rriQ and the lattice spacing a, 
so the short-distance coefficients C l£ ^ depend 
on niQa. Since all dependence on niQa is 
isolated into the coefficients, this description 



shows that heavy-quark lattice artifacts arise 
only from the mismatch of the and their 
continuum analogs Ccont- 

For methods 1 and 2, Eq. (||) is just a 
Symanzik LE£. For lattice NRQCD we re- 
cover the well-known result that some of the 
coefficients have power-law divergences.!^ So, 
to take the continuum limit one must add 
more and more terms to the action. This 
leaves a systematic error, which, in practice, 
is usually accounted for conservatively. 

Eq. ([}]) is more illuminating for meth- 
ods 3 and 4, which use Wilson fermions 
(with an improved action). Wilson fermions 
have the same degrees of freedom and heavy- 
quark symmetries as continuum QCD, so the 
HQET description is admissible for all mga. 
Method 4 matches the coefficients of Eq. (||) 
term by term, by adjusting the lattice action. 
In practice, this is possible only to finite or- 
der, so there are errors (Cj$ - Ccont )(^hqet)> 
starting with some n. Method 3 reduces mqa 
until the mismatch is of order (rriQa) 2 <C 1 
(or ^ 1). This runs the risk of reducing uiq 
until the heavy-quark expansion falls apart. 

The non-expert can get a feel for which 
methods are most appropriate by asking him- 
self what order in Aqcd/to{, is needed. For 
zeroth order, method 1 will do. For the first 
few orders, the others are needed, although 
with method 3 one should check that the cal- 
culation's Aqcd/wq is small enough too. 

3 New Results 

3.1 A and Ai 

The matching of lattice gauge theory to 
HQET provides a new way to calculate mar 
trix elements of the heavy-quark expansion Ep 
The spin-averaged B*-B mass is given bjo 

M = m + A-Ai/2m, (4) 

where m is the heavy quark mass, and M = 
|(3M B . + M B ). The lattice changes the 
short-distance definition of the quark mass:E3 

Mi -mi= Aiat - Ai lat /2m 2 . (5) 
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Because the lattice breaks Lorentz symmetry, 
mi ^ 77i2, but thexare still calculable in per- 
turbation theoryJla The lambdas in Eq. (||) 
are labeled "lat" because they suffer lattice 
artifacts from the gluons and light quark. 

After fitting a wide range of lattice data 
to En, (|^) and taking the continuum limit, we 
findtl A = 0.68±£?2 GeV and Ai = -(0.45 ± 
0.12) GeV 2 in the quenched approximation. 
The lambdas appear also in the heavy-quark 
expansion of inclusive decays. Although the 
current analysis is thorough, there are several 
ways to improve it.L2l For example, Ai at has 
an unexpectedly large a dependence, so the 
analysis should be repeated with an action 
for which C ct f in Eq. (§) is 0(a). 

3.2 B — > nlf form factors and V u b 

It is timely to discuss B° —> 7r + l~i>, be- 
cause there are three calculations to coni= 
pare, using lattice NRQCD (method_2),ll2 
the extrapolation method (method 3), til and 
the HQET matching method (method 4). 113 
UKQCD's workJsJinaljLZl and the other two 
are preliminaryOlIS The decay rate requires 
a form factor, called f+(E), which depends 
on the pion's energy in the B's rest frame, 
E = v-pir- It is related to the matrix element 
(tt\V^\B) : which can be computed in lattice 
QCD. The systematics are smallest when the 
pion's three-momentum is small. 

The three recent results are compared 
in Fig. ^. The error bars shown are sta- 
tistical only. For NRQCD these are larger 
than expected.E3 For the other two, the com- 
parison gives a fair idea of systematics that 
are not common to both, because the ex- 
trapolation of heavy quark mass needed with 
method 3 amplifies the_statistical error.113 
The other two workstStLS compute directly 
at the b quark mass and, thus, circumvent 
this problem. The heavy quark masses of 
UKQCD0 are all below 1.3 GeV, and as low 
as 500 MeV, so one might worry whether the 
heavy-quark expansion applies. 



3.3 B~ -> D 0( *h-9 and V cb 

The form factors J- B _> D (*) of the decays B — > 
D^lv are normalized to unity for infinite 
quark masses. What is needed from lattice 
QCD, therefore, is the deviation from the 
unity for physical quark masses. Hashimoto 
et al&BtB have devised methods based on 
double ratios, in which all the uncertainties 
cancel in the symmetry limit. Consequently, 
all errors scale as T — 1 , not as T . 

For B Dlv they findEI (published) 

Fb^d{1) = 1.058±0.016±0.003l^5, (6) 

where error bars are from statistics, adjusting 
the quark masses, and higher-order radiative 
corrections. For B — » D*lv they findc3 (still 
preliminary) 

•F. b—>d* (1) = 0.935 ± 0.022±°;™s 

±0.008 ±0.020, 1 ' 

where now the last uncertainty is from 1 /toq . 
In both results, an ongoing test of the lattice 
spacing dependence is not included, but that 
will probably not be noticeable. More seri- 
ously, these results are, once again, in the 
quenched approximation, but the associated 
uncertainty are still only a fraction of T — 1. 
Both results will be updated soon, with cal- 
culations at a second lattice spacing and re- 
finements in the radiative corrections. 
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Figure 3. Recent results for the decay B — > nlu. 
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group method £? Bd (4.8 GeV) 



fp( n f) 



MIL 



w. 



CP-PACSE^ 



JLQCI 

APES 
UKQCB-N 



0.85 ±0.03 ±0.11 



0.93 ±0.08 
0.92 ± 0.04l°3 



+00 
-06 



Table 1. Recent quenched results for Bg. 

3.4 Bq-B® mixing: Bb, js, and Vt q 
The mass difference of CP eigenstates is 



16tt 2 



\v;m z mq« b = 2 ), 



'tq 



Am B o 

(8) 

where the light quark q is d or s, Sq is an 
Inami-Lim function, and (Q^ B=2 ) is 

(B q \Q^\B a q )^lm% g fl q B Bq . (9) 

The \x dependence in tjb and Q^ B=2 cancels. 
New physics could compete with the W and 
t box diagrams and change Eq. (@). 

Lattice QCD gives matrix elements, so 
the basic results are {Q^ B=2 ) and /b . It 
is often stated that uncertainties in Bs q and 
£ 2 = f B BbJ f B Bb s should be small, be- 
cause they are ratios. Some cancelation 
should occur, but only if one can show that 
the errors are under control. At present there 
are unresolved issues in method 3, so one 
should be cautious. 

With that warning, results for Bb from 
three groups are in Table @. Note that 
JLQCD now includes the short-distance part 
of the 1/toq contributionE3 APEEa (final) 
and UKQClJEI (preliminary) both extrap- 
olate linearly in I/iuq from charm (e.g., 
1.75 GeV < m« B » < 2.26 GeV for APE). It is 
not clear whether the first term of the heavy- 
quark expansion is adequate here; everyone 
working in B physics can and should form his 
or her own opinion. It is also not clear how 
the (m<ga) 2 lattice artifacts of method 3 fare 
through the I/tuq extrapolation. It is likely 
that the systematic error is not well con- 
trolled and, thus, possibly underestimated in 
for the last two rows Table |l|. At present one 
should prefer the JLQCD results. 



/ B (0) 171 ± 6 ± 17+^ 

(2) 190±6± 2 °±° 

/ B ,(0) 197 ± 5 ± 23± 2 | 

(2) 218 ±5^6+11 

f D (0) 199 ± 6 ± 12 ± 14 

(2) 213±4±if±S 

f Ds (0) 222±5t} 9 7 ± 1 5 

(2) 240±4±»±g 



190±3±9 
215 ±11 ±11 
224±2±15 
250 ± 10 ± 13 
224±2±15 
236 ± 14 ± 14 
252 ± 1 ± 18 
275 ± 10 ± 17 



Table 2. Preliminary unquenched results for fg, etc., 
in the continuum limit. All values in MeV. 



The MILCEI and CP-PACSO groups 
have new, preliminary unquenched calcula- 
tions of the heavy- light decay constants fs, 
fs s , Id, and frj s . Both use method 4. Both 
have results at several lattice spacings, so 
they can study the continuum limit. The sta- 
tus for Osaka is tabulated in Table ||. The 
first error is statistical, the second system- 
atic. MILC also provides an estimate of the 
error from quenching. (With n/ = 2 the 
strange quark is still quenched.) CP-PACSc3 
also has results with method 2, which agree 
very well with method 4. One should not, 
at this time, take the differences between the 
two groups' central values very seriously. It 
is more important to understand the different 
systematics of methods 2, 3, and 4. 

4 Prospects 

For B physics it is important to remove the 
quenched approximation, more so than to re- 
duce the lattice spacing much further. To 
do so, we need more computing. Fermilab, 
MILC, and Cornell are building a cluster of 
PCs to tackle the problcm.il Our pilot clus- 
ter has 8 nodes with a Myrinet switch. We 
plan to go up to 48-64 nodes, and then hope 
to assemble a cluster of thousands of nodes. 
The large cluster would evolve, by upgrading 
a third or so of the nodes every year. This 
is an ambitious plan, but not more ambitious 
than the experimental effort to understand 



5 



flavor mixing and CP violation. 

The last few years have seen significant 
strides in understanding heavy quarks in lat- 
tice QCD. The progress has been both com- 
putational and theoretical, with one guiding 
the other. Calculations shown here, for B a — > 



7T+Z-P, D° 



'l + v, B~ -» D°^l-9, and 



Bj-Bj mixing are a subset, but in Fig. gthey 
are as basic as A, B, C. With the right 
amount of support from the rest of the com- 
munity, we hope to obtain the tools needed 
to resolve the few outstanding problems and 
to produce excellent unquenched results. In- 
deed, the example of fs shows that this is 
already beginning. 
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Figure 4. Artist's conception of the author striving to comprehend charm and beauty. (©2000 Mercedes 
Kronfeld Jordan.) 
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